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TECHNICAL MEMORANDUM 


COMPARISON OF TWO COMPUTER CODES FOR CRACK GROWTH ANALYSIS — 

NASCRAC VERSUS NASA/FLAGRO 


The structural integrity of space flight hardware is established by a combination of qualifica- 
tion tests and analyses which simulate actual operating conditions, including flight loads, tempera- 
tures, and corrosive environments. These structural analysis and test activities usually fall into three 
distinct areas. The first two areas, strength and fatigue analysis, assume the load carrying structure 
is unflawed. This assumption implies that no defects have been introduced during the manufactur- 
ing process of each individual part, which in reality can never be possible on an economical basis. 

The existence of flaws is accounted for in the third area, fracture mechanics. This aiea 
becomes an important effort in which defects are known as a result of quality inspections, or 
assumed to exist in a part and an assessment is made as to their impact on the parts useful life. 
Fracture mechanics attempts to predict the useful service life of an initially flawed structural part 
by calculating crack growth and eventual part failure due to unstable crack growth. 

This paper compares the service life calculations of two computer codes, NASCRAC and 
NASA/FLAGRO. The analysis technique is based on linear elastic fracture mechanics (LEFM), in 
which stresses remain below the yield strength of an elastic/plastic material. Subcritical crack 
growth calculations assume that in a metallic part, the extent of yielding at the crack tip is very 
small compared to the crack size, uncracked ligament, and the bulk of the cracked body remains 

elastic. 


To perform service life calculations, one must have a relationship expressing incremental 
crack growth, DA/DN, as a function of loading, geometry, and crack size. Load, crack size, and 
geometry are expressed in terms of the cyclic stress intensity factor, AK. The crack growth iate as 
a function of AK is then determined by material tests, plotting DA/DN versus AK for the given 
material form and H.T. condition, loading condition, and environment. 

Crack growth rate equations such as the Paris, Walker, and modified Forman equations are 
used to obtain a “best fit” curve to the laboratory DA/DN versus AK data. Constants in the 
equations which result in a “best fit” then become crack growth rate material constants for a 
particular set of laboratory conditions. 

Two extreme values of AK also become material constants; AK 0 is the threshold stress 
intensity below which no crack growth occurs, K c is the critical stress intensity at which a crack 
becomes unstable and complete fracture occurs. Formulations of AK solutions and crack growth 
rate equations form the basis of computer codes which numerically integrated the 
DA/DN = F(AK) relationship. 

Before a computer code is used as part of the structural integrity assessment process, it 
should be exercised thoroughly and its numerical calculations checked to insure reasonable and 
accurate answers. The results presented herein compare the Safe Life calculations of two computer 
codes with each other, and with test data to a limited extent. 



The computer program NASA/FLAGRO (commonly known as NASGRO) became available 
in 1986 from the NASA Johnson Space Center. The program was developed under the guidance of 
the NASA Fracture Control Analytical Methodology Panel and contains stress intensity factor 
solutions to a number of commonly used crack geometries. Service life calculations are performed 
with the modified Forman equation which reduces to the Walker or Paris equation depending on 
material constants used. 

NASA/FLAGRO is menu driven and prompts the user for information in a serial manner. 
After selecting the type of analysis desired, such as Safe Life, the user answers a series of ques- 
tions and enters data depending on the particular path taken. Generally, the program operates 
serially, requiring the user to follow the same path and answer a number of basic questions before 
each execution. 

The computer program NASCRAC is being developed by Failure Analysis Associates under 
contract to NASA at Marshall Space Flight Center. For Safe Life analysis, NASCRAC has basically 
the same capabilities as NASA/FLAGRO, although implemented differently. Generally, stress 
intensity factors are obtained from influence function solutions to various geometries for which 
exact solutions do not exist. NASCRAC enables the user to select any one of several commonly 
used crack growth equations; including the Paris, Walker, and modified Forman equations. 

NASCRAC is similar to NASA/FLAGRO in that the program is menu driven and the user 
answers questions and enters data in response to screen prompts. With NASCRAC, however, the 
user is not required to answer a series of questions before each execution. The user may randomly 
select only those menu items relating to the particular solution desired. 


NASA/FLAGRO SAFE LIFE FEATURES 


NASA/FLAGRO features which can affect Safe Life calculation: 

(1) For surface cracks with constant amplitude loading, AK is multiplied by a crack closure factor 
Pr. 


Pr — 


0.9 + 0.2 R 2 - 0.1 R 4 
0.9 


; R > 0 
; R 0 


This can increase fatigue life. 

(2) AK lh , the fatigue threshold is calculated using, 


AK, h = (1 - C 0 R) d AK U 
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To be conservative, let C G — d = 1 for R 2= 0 


AK th - (1-R) AK 0 

For small cracks, a ^ 0.025 in, AK 0 =• 0. 

(3) Input of 

K IC - plane strain fracture toughness 

K Ie - fracture toughness for an elliptical crack 

A k , B k - fit parameters 

To calculate K c - critical stress intensity, 



c) K c — K IC (1 + B k e- w ) 

K c is incorporated into the modified Forman equation to accelerate DA/DN as K c is approached. 

NASCRAC SAFE LIFE FEATURES 

(1) Crack growth equations 

a) Modified Forman Analytical Comparison 

b) Walker Comparison to test data on both codes 

(2) Piecewise Linear Approximation method used. 

(3) K-solutions are based on influence functions with the default order of accuracy. 



ANALYSIS NOTES 


Surface Flaws 

NASCRAC — Uses K lc value to accelerate DA/DN per the Forman equation and 
defines failure when AK > K] C , where K| C is manually input. 

NASA/FLAGRO — Uses K c value calculated from K) C and other variables to 
accelerate DA/DN per the Forman equation when AK > K le , where K] e is a material constant 
for surface flaws. 


Growth Rate Equations 

NASCRAC — Uses the following growth equations: Paris, Modified Forman, Walker, 
Collipriest, and Hop Rau. 

NASA/FLAGRO — Primarily uses the Modified Forman Equation but the Paris and 
Walker equations could be used. 

K c Values 

NASCRAC — K l . is used in the Modified Forman equation but K is the controlling cutoff 

value. 

NASA/FLAGRO — For B k =f 0, NASA/FLAGRO uses a K lc larger than K, c for thin 
material (Reference 8) 

K,./K lc = 1 - B k e- W 

when B k = 0, K c = K lc , where w = (A k t/t„) 2 . 
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COMPARISON ANALYSIS CHART 


Type of Geometry 


Through Center Crack 


Through Edge Crack 


*See analysis results. 


Parameters 

Type of Run 

NASGRO 

NASCRAC 

W = 10.0 

t = 0.25 

*R = 0 Tension Only 

X 

X 

4130 Steel 

a ; = 0.05 

*R = -1 Tension Only 
Closure 

X 

X 

or t - 50 Ksi 

No Closure 

X 

X 

W = 10.0 
t - 0.25 

4130 Steel 

*R = 0 Tension Only 

X 

X 

cr, = 50 Ksi 

o b — 50 Ksi 

*R = -1 Tension Only 
Closure 

X 

X 

a j = 0.05 

No Closure 

X 

X 


*R = 0 Bending Only 

X 

X 


*R = — 1 Bending Only 
Closure 

X 

X 


No Closure 

X 

X 


*R = +0.5 Tension 
Closure 

X 

X 


No Closure 

X 

X 


*R = - 0.5 Tension 
Closure 

X 

X 


No Closure 

X 

X 


R = +0.5 Bending 
Closure 

X 

X 


No Closure 

X 

X 


R = -0.5 Bending 
Closure 

X 

X 


No Closure 

X 

X 

a = 0.25 

R = 0 Tension 

X 

X 


R = -1 Tension 
Closure 

X 

X 


No Closure 

X 

X 
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Type of Geometry 

Parameters 

Type of Run 

NASGRO 

Through Crack at 
Pin Loaded Hole 

W - 1.75 
t = 0.44 
D = 0.375 
B - 0.83 
4340 Steel 
cr t = 59 Ksi 

*R = 0 Tension 

X 

Through Crack at 
Pin Loaded Lug 

o"b — 37 Ksi 
a - 0.25 

W = 5.0 
t = 0.25 
D = 0.5 
a t = 150 Ksi 
4130 Steel 
a| = 0.05 

+ 

Bearing 

R — 0 

X 


a; = 0.10 

R = 0 

X 


a ; = 0.25 

*R — 0 

X 

Surface Flaw Center 

Crack Specimen 

Test Spec. No. 

62* 



W = 4 
t - 0.50 
cr t = 84 Ksi 

aj = 0.06 
aj/2Cj = 1/2 

Ti = 6AL-4V 

R - +0.05 

X 


Test Spec. No. 
W = 4 
t = 0.50 
CT t — 57 Ksi 

aj = 0.06 

a i /2c i = 1/2 

5576* 

R = +0.05 

X 


Ph- 1 3-8M c 




*See analysis results. 


NASCRAC 


X 


X 
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XXX 



NA3QRQ 



TC02 


GEOMETRY 
MODELS HASCRAC 



THROUGH 

CENTER 

CRACK 



202 


THROUGH 

EDGE 

CRACK 



203 


Sq ♦P/Wt 

t t t t f t 


THROUGH CRACK 
AT HOLE 
PIN LOADED 

So 

Si «P/Dt 
TC03 



»»!»»» 


zz 



°zz 

208 
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MASQRQ 


HAScrac 


p/wi 



THROUGH CRACK 209 

AT LUG 



SC01 


Z 



SURFACE CRACK 
(CENTER) 


y 



702 


B 




ANALYSIS RESULTS 
THROUGH CENTER 
CRACK 

4 w-io” 

1 t==0 - 25 ” 

“I 2 aj = 0.05” 

4130 STEEL 

(7 = 50KSI 

1) R = 0 

2) R = 1 

i 



NASGRO MODEL TYPE TC01 NASCRAC MODEL TYPE 202 


FIRST CASE R=Q 


NASGRO 

k MAX = 80 - 10 KSI'Hn @ 20,173 CYCLES 2a=1.584" 
NASCRAC 

K max =80 KSifTN @ 20,176 CYCLES 2a=1.60" 


SECOND CASE R=-1 
NASGRO 

CLOSURE: K MAX =80.12 KSI flN @ 16,401 CYCLES 2a=1.584" 
NO CLOSURE: K MAJ f=80.04 KSI fTFl @ 4,459 CYCLES 2a=1.58" 

NASCRAC 

CLOSURE: K MAX =80 KSHn @ 4433 CYCLES 2a=1.60" 

NO CLOSURE: SAME AS CLOSURE 
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THROUGH EDGE 
CRACK 


Si 



t So 


a 


siC^ 


I So 

■^3 


W s 10" 
t = 0.25" 

C| = 0.05 
4130 STEEL 
0 t = 5OKSI OR 

1 ) R = 0 

2) R = -1 

3) R = +0.5 

4) R = -0.5 



a b = 50 KSI 

DRtO 
2) R s-1 


NASGRO MODEL TYPE TC02 NASCRAC MODEL TYPE 203 

TENStON_.ONLY 


FIRST CASE R = 0 
NASGRO 

K MAX s 80 KSI ® 9674 CYCLES a f= °* 610 " 


NASCRAC 

K max = 83.22 KSI JW <§> 9616 CYCLES a f = 0.655" 


SECOND CASE R = -1 



NASGRO 

CLOSURE: 

k MAX = 80 KSI @ 7901 CYCLES 

a f * 0.61” 

NO CLOSURE: 

k MAX 58 80 KSI ® 2148 CYCLES 

a f s 0.61" 

NASCRAC 

CLOSURE: 

k MAX = 83 20 KSI @ 4439 CYCLES 

a? = 0.655 

NO CLOSURE: 

SAME AS CLOSURE 
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THROUGH EDGE CRACK 
TENSION CASES CONTINUED 


THIRD CASE R=+0.5 


NASGRO 


CLOSURE K max = 80.05 KSlV IN @ 27,485 CYCLES a=0.61" 


vuvwwi li- ''MAX — ” - 

NO CLOSURE K max = 80.05 KSlY IN @ 43,772 CYCLES a=0.61 


NASCRAC 

CLOSURE K max =83.22 KS1 1 1N @ 43,516 CYCLES a=0.655" 
NO CLOSURE SAME AS CLOSURE 


FOURTH CASE R=-0.5 
NASGRO 

CLOSURE K max = 80.00 KSI YIN @ 8 757 CYCLES a=0.61" 
NO CLOSURE K max = 80.00 KSI /Tfr@ 4010 CYCLES a=0.61 


NASCRAC 

CLOSURE K max = 83.22 KSMMN @ 3985 a=0.655" 
NO CLOSURE SAME AS CLOSURE 



THROUGH EDGE CRACK CONTINUED 
BENDING ONLY 


FIRST CASE R=Q 
NASGRO 

K MAX = 80 - 09 KsiVTrT @ 10,830 CYCLES a =0.73 in. 

NASCRAC 

k max = 81 - 363 KsiVTiT@ 10,228 CYCLES a =0.721 in. 

SECOND CASE R=-1 
NASGRO 

CLOSURE: K MAX =80.03 KsiVTrT @ 8846 CYCLES a =0.73 in. 

NO CLOSURE: K MAX =80.03 KsiVTrT @ 2405 CYCLES a =0.73 in. 

NASCRAC 

CLOSURE: K MAX =81.36 KsiVTrT @ 2271 CYCLES a =.721 in. 

NO CLOSURE: SAME AS CLOSURE 
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THROUGH CRACK 
AT PIN LOADED HOLE 


W 



W a 1.75" 
t = 0.44" 

HOLE DIAMETER a 0.375" 
EDGE DISTANCE a 0.83” 
4340 STEEL 

(7 T a59 KSI + 0 bear =37 KSI 
CRACK LENGTH a 0.05" 

R a 0 


NASGRO MODEL TYPE TC03 NASCRAC MODEL TYPE 208 
NASGRO 

K MAX = 90 * 17 KSI ® 4 ’ 334 CYCLES a f =0-339" 
NASCRAC 

k MAX = 90.13 KSI 'Hn <§> 6609 CYCLES a f =0.492" 




THROUGH CRACK 
AT LUG 



P/wt 


NASGRO MODEL TYPE TC04 NASCRAC MODEL TYPE 209 

WIDTH = 5.0” 

THICKNESS = 0.25” 

4130 STEEL 

DIAMETER OF HOLE 0.5” 
a, = 0.25” 

<r = 150KSI 


NASGRO RESULTS 

K max = 80.71 KSlflNf @ 64,426 CYCLES a/il.99” 
NASCRAC RESULTS 

K M AX=80 KSI iffFT" @113,649 CYCLES a/= 2.184" 



WALKER CONSTANTS FOR PART 
THROUGH CENTER CRACK ANALYSIS 


MATERIAL: PH 13 - 8 M 0 TEST CASE N0 - 5576 


c 

7.63 x 10’ 11 IN/CYCLE 

m 

1.0 

n 

3.54 

AK th 

8 KSI V IN 

K,c 

100 KSI V IN 

MATERIAL: 

Tl - 6AL - 4V TEST CASE NO 

c 

2.914 x 10’ 12 IN/CYCLE 

m 

0.04435 

n 

4.51 

AK t (, 

4.5 KSI Y IN 

K .c 

70 KSI V IN 
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a or c 


NASA/FLAGRO PART THROUGH CRACK ANALYSIS 


DATA FROM "NASGRO #5576" 



TRANSITION TO 1-D SOLUTION TC01 AT 43,531 CYCLES 

a=0.5" t=0.5" c=0.62" 

FAILURE OCCURED AT 43,843 CYCLES K max =100.3 KSlIfiN c=0.798” 








DEP TH a QR H ALF L ENG T H C 



CO 



NASA/FLAGRO PART THROUGH CRACK ANALYSIS 

DATA FROM "NASGRO #62" 



2000 


4000 


6000 


CYCLES 

FAILURE OCCURED AT 7855 CYCLES a=0.358" C=0.4" Kmax=85.09 KSI 



NASCRAC PART THROUGH CRACK ANALYSIS 


DATA FROM "NASCRAC TEST 62" 


H-2CX 


T 


WIDTH=4 IN. 
tsO.5 a=0.06 a/2c=1/2 


CYCLES 


FAILURE OCCURRED AT 7214 CYCLES 
Kmax - 71 KSIVlN = 0.367” 




OBSERVATIONS 


NASCRAC and NASA/FLAGRO are both user-friendly fracture mechanics analysis codes. 
Both programs offer a wide variety of crack geometries. Material property data can be read in from 
a resident file or from user defined input. Load spectra data for the constant amplitude loading 
cases were utilized easily in both programs. 

For the through-crack comparison analysis the Modified Forman equation was used and for 
the part- through crack analysis the Walker growth rate equation was used. 

For the through-crack analysis with an R ratio of zero, results showed good correlation 
between the two codes except for the through-crack at a lug solution. For R = -1, +0.5, -0.5, 
NASA/FLAGRO calculates an m value that is not readily known to the user; it must be hand 
calculated for use in NASCRAC. By specifying the nonclosure option, m is automatically set to 
zero. The nonclosure option gave the most conservative results in NASA/FLAGRO. For R = -1, 

+ 0.5, —0.5, changing the m value in NASCRAC had no effect on the results, the m value has been 
permanently set to some prescribed value. The NASCRAC results for the through-crack analysis for 
R = -1, +0.5, -0.5 were in the range of the NASA/FLAGRO results for the nonclosure option. 

For the part through center crack analysis, both programs gave comparable results, particu- 
larly with specimen No. 15576 where the crack grew through before failing, but both programs 
showed failure before breakthrough for specimen No. 62 which was different from the results of 
the test. 

The comparison analysis between the two programs is an on-going effort for our analysis 
team. Other types of solution methods and problems are scheduled to be studied in the future. 
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